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Abstract 
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matrix sampled from the real Ginibre ensemble has 2k real eigenvalues. We 
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1 Introduction and the main result. 


Our paper is dedicated to the study of the probability p 2 n, 2 k that a real 2n x 
2 n random matrix with independent normal entries (the so called ‘real Ginibre 
matrix’) has 2k real eigenvalues. It has been known since |8j that a typical large 
N x N Ginibre matrix has 0{y N) real eigenvalues. What is the probability of 
rare events consisting of such a matrix having either anomalously many or few real 
eigenvalues? 

The former question has been addressed by many authors. Building on the 
original work [12] by Ginibre, Edelman used the real Schur decomposition to prove 
that 


Pn,n = 



N(N-l) 

4 


see PJ. In [5j Akemann and Kanzieper employed the method of skew-orthogonal 
polynomials to determine the probability that all but two eigenvalues of a real 
Ginibre matrix are real. In the large -N limit, their result reads 

„ _ ^21 N 2 + ^sip^l N+o{N) m 

Pn,n-2 — e 4 2 - , fij 


where lini/v->.oo o(N)/N = 0. These answers were generalized in a very recent pa¬ 
per [TS] where the large deviations principle of [3] was extended to prove that the 
probability that a real Ginibre matrix has aN (where 0 < a < 1) real eigenvalues 

is PN,aN A ~°° e -Ar2/ “, where the symbol denotes the logarithmic asymptotic 
equivalence and the constant I a is characterised as the minimal value of an explic¬ 
itly given rate functional, see Proposition 2 and formula (4) of [13]. 

In the present paper we answer the question about the probability that a real 
Ginibre matrix has very few real eigenvalues: 

Theorem 1. Let G 2 n be a random 2 n x 2 n real matrix with independent N(0, 1) 
matrix elements. Let p 2 n, 2 k be the probability that G^n has 2k real eigenvalues. 
Then for any fixed k = 0,1, 2, 3,..., 

lim 4= logP2n,2fc =-4=C ( | ) , ( 2 ) 

V2n V27T \2 J 

where ( is the Riemann zeta-function. Moreover, 

h m »^ logP2 "* = -^ c (^)’ (3) 

where (fc n ) n >i is a sequence of non-negative integers such that 

lim (n -1 / 2 log(n)) k n = 0. 

n —>oo v 7 
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In particular, the probability that a large 2 n x 2 n Ginibre matrix has no real 
eigenvalue behaves as 


p2n,0 e -V!W»^>. 

Notice that the answer (|2]) is qualitatively different from the results for the 
probability of having 0(n ) real eigenvalues quoted above: the ‘cost’ of having 
0(n ) real eigenvalues normalised by the total number of ‘anomalous’ eigenvalues 
increases linearly with n, whereas the ‘cost’ of removing all real eigenvalues from 
the real axis is constant per eigenvalue. 

It is also worth noting that our result ‘almost’ extends to the typical region 
k ~ n 1/2 . (For example, we can choose k n = \\fn/ log 2 n] in ^.) It would be 
interesting to see if (J3]) survives for k n = \csjn\ where c « 1. 

The statement of the theorem can be guessed using existing results: in the limit 
N —> oo, the unsealed law of real eigenvalues for the real Ginibre N x N ensemble 
converges. The limit coincides with the t = 1 law for the A + A —> 0 interacting 
particle system on M [16]. The probability that an interval of length s has no 
particles for A + A —> 0 has been calculated formally by Derrida and Zeitak [7], 
These two facts allowed Forrester m to conclude that the large-A" limit of the 
probability that there are no real eigenvalues in the interval (a, a + s) should be 
given by 


Prob[Goo has no eigenvalues in (a, a + s)] ~ e 2 V 2 V^ 3//2 ^. ( 4 ) 

Let us stress that equation Q is valid for N = 00 only. However, we know from the 
work of Borodin and Sinclair [ 6 j and Forrester and Nagao [TO] that the law of real 
eigenvalues for the real Ginibre ensemble is a Pfaffian point process for all values of 
N ^ cc. Convergence of the finite -N kernel to the N = co kernel is exponentially 
fast within the spectral radius. The spectral radius is Rn = VN + 0(1) [ 8 ]. We 

also know that the boundary effects for a large but finite matrix size N are only 

felt in the boundary layer of the width of order 1 near the edge. Therefore, the 
simplest finite -N guess for Prob\GN has no real eigenvalues] is 

Prob\G]y has no real eigenvalues] 

% Prob[GN has no real eigenvalues in (—Rn + L,Rn — L)] 

% Prob[G x has no real eigenvalues in (— R N , Rn)] 

Here L » 1 is a large N -independent constant. The last probability in our heuristic 
chain of arguments can be approximated using Q with s = 2 Rn■ This suggests 

Prob[GN has no real eigenvalues] % , 
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Figure 1: The logarithm of the probability p ^,o that an N x N matrix of even size 
sampled from the real Ginibre ensemble does not have any real eigenvalues, as a 
function of yN. The leading coefficient extracted using the best fit is —1.042, the 
best fit for the next-to-leading constant is 0.06268. The ‘exact’ curve is constructed 
using formula (|6]) of Lemma [I] below. The form of the & 2 -term in the fitting curve 
was chosen to minimize the numerical goodness-of-fit y 2 . 


which agrees with the statement of Theorem [l] 

The value of the constant which defines the rate of decay of p 2n> o in (J2| is 


1 


V2vr 


C(3/2) % 1.0422, 


which is consistent with its numerical estimate, see Figure [l] The numerical analy¬ 
sis of the exact formula for p 2n ,o (see (|6]) below) also shows that under the assump¬ 
tion that the next-to-leading term in the large-lV expansion of p ^,o is constant, 
the resulting coefficient (% 0.06268) is close to its exact counterpart from the large 
gap size expansion of the Derrida-Zeitak formula (% 0.0627). At the moment we 
do not have a theory explaining this closedness. 

Both the numerical simulations and the heuristic argument given above provide 
a strong hint in favour of Theorem [lj 

There are several possible routes to the proof of the theorem. For example, one 
can try to use Forrester’s observation, coupled with the knowledge of the rate of 
convergence of the Borodin-Sinclair-Forrester-Nagao kernel in the large -N limit, to 
show that the errors in applying Derrida-Zeitak’s formula to gaps of iV-dependent 
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sizes vanish as N —* oo. There is however a problem with this approach: in 
the case we are interested in (annihilating Brownian motions or the 2-state Potts 
model) the infinite sums entering the gap formula converge only polynomially, 
see [7] for details. Therefore, a careful justification would be required for the 
validity of the interchange of summation and taking the large gap size limit. We 
feel that such a justification is best done in the context of a general theory of 
‘Fredholm Pfafhans’. In this paper, we will adopt the spirit of Derrida-Zeitak’s 
calculation to construct rigorous asymptotics of a very compact and easy to use 
exact determinantal expression for the probability P 2 n , 2 k specific to the real Ginibre 
ensemble. This determinantal expression can be derived building upon the results 
of [13] and [IB], see Lemma [l] below. We hope of course that our very specialized 
proof will contribute to the general discussion of the theory of large deviations for 
Pfaffian point processes. 

There is a drawback to our approach as well: even though we can now claim 
that (]2]) is true, we still do not know how a large Ginibre matrix without real 
eigenvalues looks. For example, is there a unique optimal configuration of complex 
eigenvalues for such matrices? What can be said about the overlaps between left 
and right eigenvectors of Ginibre matrices without real eigenvalues? To answer 
these questions, one has to develop a large deviations principle along the lines of 
na which will most likely use the picture of the ‘two-component’ plasma consisting 
of one-dimensional and two-dimensional ‘gases’ of eigenvalues discussed there. 

Our paper is organised as follows: a reader who is satisfied by our heuristic 
argument and the numerics can stop here. Those interested in the mathematical 
proof are advised to read Section 2 and consult Appendix [A] for the proofs of the 
technical facts used in the proof of Theorem [T[ Appendix [B] contains remarks on 
the numerical evaluation of p 2 n,o f° r large values of n. 


2 The proof of Theorem [T|. 


Our starting point is the following exact determinantal representation for the gen¬ 
erating function for the probabilities p 2 n, 2 fc : 


Lemma 1. Let n be a positive integer. Then 

{z ~ 1) r (j + k- 3/2) 


k=0 

In particular, 


y z k p 2n , 2 k = det 

j,k=l,n 


8j,k + 


V2/r y/T(2j - l)T(2/c - 1) 


P2n,o = det 

j,k=l,n 


3j,k 


F(j + k - 3/2) 


V 2 tt - l)r(2fc- 1) 


(5) 


( 6 ) 
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We postpone to Appendix A the proofs of all lemmas used during the proof of 
the main theorem. 

Notice that the expression (| 6 ]) coincides (as it should) with the s —* go limit of 
the probability that a 2 n x 2 n real Ginibre matrix has no real eigenvalues in the 
interval (— s,s) calculated by Forrester, see formula (3.58) of [IT ]. 

We will prove Theorem |T] in two steps: first, we will prove (j2j) for k = 0, then 
we will show that lim^* 7 = logp 2 n, 2 fc„ = lim^* -R logp 2 n,o, where (k n ) n>1 is a 
sequence of integers which grows with n slower than n 1 ' 2 / log(n). 


2.1 The calculation of lim n ^oo logp 2 n,o* 

Let M n be an n x n symmetric matrix entering the statement of Lemma [l] 


M n (j, k ) 


i r (j + k - 3 / 2 ) 1 , 

—, 1 < j, k n 

V2^V r ( 2 J - WW- 1) 


(7) 


Lemma 2. M n is a positive definite matrix. Moreover, there exists a positive 
constant p > 0 and a natural number N such that for any n > N, 


n 


where A max (n) is the maximal eigenvalue of M n . 

Using Lemmas [T| and [ 2 ] we represent p 2 n,o as follows: 

1 logp 2 n,o = ~ 7 =Tr log(J — M n 


V 2 


n 


V2n 


K n 


V _ TrAC - - 7 = 
^/2n “ m ^fi2n 

v m =1 v 


Rn(K n ), 


( 8 ) 


(9) 


where K n is a cut-off which increases with n (chosen below) and R n is the remainder 
of the Taylor series for log(J — M n ) written in the integral form: 

r 1 / M K+l \ 

*■<*>-/. ^( d 

An upper bound on \R n (K)\ follows from Lemma [2] by replacing all eigenvalues of 


6 










M n with A max(n): 


\Rn{K)\ ^ n\«: x \n) 

< nX maxH l ° g 


(1 ~x) 


K 


0 (1 - A max{n)x) K+l 

1 


dx 


1 — A 

< n\og ( - I (1 - 


n 




max (n) 
K 


TlJ 


So, if we choose 


it is easy to check that 


K n = [n a \, a > 1 , 


lim R n (K n ) = 0. 


( 10 ) 


( 11 ) 


The last step of the proof is the calculation of X! m =i “TrM™. The relevant 
results can be summarized as follows: 


Lemma 3. For any fixed integer rn > 0, 


lim _ 

n ->oo V2 n 


Tr M™ = 


277171 


( 12 ) 


Moreover, for any any positive integers m, n 


m 


TrM™ < ./— (1 + n- 1 ) + 7 + h/— (1 + 2rr‘) 
x/ 7rm v y 4 8 V ra v ' 


(13) 


Let us stress that formula (12) alone is not enough for the calculation of the 
lim^oo n~R 2 logp 2 n,o using § since the limits n —> oo and m —> oo do not neces¬ 
sarily commute. Instead, let us fix an arbitrary integer K > 0. For a sufficiently 


large n (so that K n > K) relation ([13) gives: 

1 Kn 1 

2 -TrM^ 


l ^ l 


-i= Y! — TrM n ^ / J 

\/2n i m \[2m m 

v m=l v m =1 

(1 + n" 1 ) - 


V 2 vr 


^ m " 3 / 2 + 

m=l 


l y i 


2 


-b 


4\/2n “[ m 8\/27r 


,-r* 

2 


(1 + 2 n~ l ) 1 


71 “i' /m 


(14) 


In writing the above double inequality we used the fact that M n is positive definite, 
which implies that Tr M™ > 0 for all values of m, n. Let us choose K n in the form 
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(10) with a < 2 and take n —* oo in (14). As K is n-independent, we can use 
formula (12) to compute the limit of the left hand side. On the right hand side, 
the last two sums vanish in the limit (as log [n)/^/n and n a ^ 2 ~ l correspondingly). 
The first sum converges to 


— y 

V2^r 


m 


-3/2 _ 


■\f2n 


C(3/2), 


where ((x) = X!m=i m x is the Riemann zeta-function. 
We have found that for any positive integer K , 


1 K 

— y 

\Z2tt ^ J 1 

v m= 1 


m 


-3/2 


i y i 


^ lim inf _ V —TrM, 

n—*co V2 n , m 

v m—1 

l l 

^ lim sup ._ V —TrM, 

n->«> V 2 n ^ m 


m 

n 




V2vr 


C(3/2). 


As A' is arbitrary, we conclude that 


^ K n 

lim _ 7 

n ~> oo V2n , 

v m= 1 


S TrM; 


VS 


C(3/2). 


(15) 


So we proved that both (11) and (15) hold provided the cut-off is taken in the 
form (10) for any fixed a e (1, 2). 

Finally, we can take the n — » oo limit in ([9]). Employing ( fTT| and (15) we find 
that 

hm ^=logp 2ni0 = —i=C(3/2). 

«—► co V2n V27T 

Theorem[l]is proved for k = 0. 

2.2 The calculation of lim n ^oo -^== log A 2 n, 2 /c for k > 0. 


It follows from Lemma Q] that 


Equivalently, 


P'2n,2k = [j:) det(J + (z - l)M n ) | z=0 • 


P2n,2k — 


P2n,0 { d 


k\ V dz 


— det (/ + zP, 


n) \z= 0) 


(16) 





















where P n = (I — M n ) 1 M n . Recall that 


det(J + zP n ) = z k e k (u), 

k =0 

where u = {yi, v 2 , ■ ■ ■, r'n) are the eigenvalues of P n and e k is the degree-A; elemen¬ 
tary symmetric polynomial in n variables [QL] , 

ekiy) = Yj v h v i2 .. .v ik . 


Therefore, 


P2n,2k = P2n,0 e fc (l/), for k = 0, 1, . . . , 71. 
Let us enumerate the eigenvalues of M n and P n as follows: 


(17) 


Ai ^ A 2 ^ ^ A n > 0, 

z'l ^ v 2 > ■ ■ ■ > v n > 0. 

By the definition of P n , zy = Note that u l is a monotonically increasing 

function of A*. Combining this remark with the spectral bound of Lemma [2j we 
get the following bound on the elementary symmetric polynomials: 


< z/ 1 fc e fc (l,l,...,l) ^ 


Ai 


1 - A, 


n K < [ - I n\ 


(18) 


Substituting (18) into (17) we obtain the following upper bound on logp 2n2fc : 


l°gP2n,2fc « lt>gP2n,Cl + \ ' 


(19) 


Next we derive a lower bound on logp 2rij2 fc. By positive definiteness, u t ^ \ and 
therefore e k (y) > e&(A). Let us £x a positive integer k. Due to (12), for any e > 0 
there is a positive integer N e such that for any n > N e 


On the other hand, 


/ Tl 

J 1 - e) ^ Tr M n ^ 



( 20 ) 


TrM„. — (Ai + ... + A fc _i) + (Afc + ... + A n ) ^ (A; — 1) + (n — A: + 1)A^, (21) 


9 







where the inequality is due to (J8| and the chosen ordering of A’s. 

Combining (20) and (21) we obtain the following bound on the fc-th largest 
eigenvalue of M n \ 


A k A’ 


e) — k + 1 


n — k + 1 


( 22 ) 


which holds for n > N e . Inequality (22) leads to the desired bound for efc(zz): 


efc(z') A= efc(A) ^ AiA 2 ... A& ^ A^ ^ 


' Vf(i- £ )-fc + i' 
n — k + 1 


Substituting this result into (17) we find that 


logp 2 n,2fc > logp 2 n,0 + k log 


-(1 — e) — k + 1 


n — k + 1 


(23) 


Combining (19) and (23) we find that 


1 1 
lim —=logp 2 n,2fc = hm — =logp 2 n,0- 
n— *00 V2n V2n 


(24) 


Relations (24) and (16) imply that formula ([2]) of Theorem [I] is proved for any 
fixed integer k > 0. 


Moreover, it is evident from (19) and (23) that the equality ( |24j ) generalizes to 

1 


lim 

n—>o 0 


/— logp 2 n,2fc n = lim —=logp 2 n,0- 
V2 n n-*co sJ2n 


(25) 


where {k n ) n ^\ is a sequence of natural numbers such that 

lim (n -1 / 2 log(n)) k n = 0. 

n —>oo v ' 

This proves the last claim of Theorem [lj 


Remark. Our proof of the k > 0 part of the Theorem is a simple consequence 
of positive-definiteness of M n , the spectral bound and the fact that Tr(M n ) n ~° 
■\fnpK. It is interesting that the proof does not rely on any detailed knowledge of 
the spectrum of M n . 
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A Proofs for the lemmas. 

A.l Lemma 1. 

To prove the lemma, we start with the exact formula due to Kanzieper and Ake- 
mann ra which expresses the probabilities P 2 n, 2 k hi terms of elementary symmetric 
functions: 


P2n,2k P2n,2n &n—k{tli ■ ■ ■ jtn—kji 


( 26 ) 


where t/s are given by 


tj = jTi 

Here A and B are 2 n x 2 n antisymmetric matrices whose entries 

-A'fc (flj—h Qk— i)r j 

Bjk (Qj—ii Qk—i)c j 

are defined in terms of skew products 


(27) 


(28) 

(29) 


(f,9) r=, / dxdye (l,+r)/2 sgn('i i-x)f(x)g(y) 

4 ./TO2 


(30) 


and 


(f,9)c = i 


d 2 z e~ {z2+ ~ z2)/2 erfc 


'Imz>0 


Z — Z 

iV2 


[f{z)g{z)-g{z)f{z)]. (31) 


Let us stress that (26) is valid for an arbitrary choice of monic polynomials qj(x) 
of degree j, provided matrix A is invertible. 


Substituting Eqs. (26) and (27) into the generating function 


ff2n(z) = 2 j zk P2n,2k 


(32) 


k =0 
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and making use of the summation formula [ 13 ] 


00 / 00 

= exp I 

e=o \j=i 

we obtain the Pfaffian representation mmu- 

g2n{z) = P 2n,2n Pf(-A~ 1 ) Pi{ Z A + B ), (34) 



see remark 1.3 of [5] justifying the transition from square roots of determinants to 
Pfaffians. Since g 2 niX) = 1, P 2 n, 2 n = (Pf(—A -1 ) Pf (zA + B)) 1 and (34) simplifies 
to 


Pf (zA + B) 
~ Pf (A + B ) ' 


(35) 


Next we will use the fact that expression (35) for the generating function does 
not depend on a particular choice of monic polynomials qj(x) in (28) and (29) 
to simplify it even further. Namely, we will choose qj(x) , s in such a way that 
the matrix A + B is block diagonal. Clearly, such polynomials should be skew- 
orthogonal with respect to the skew product 


(f,9) = </.J)r+ </.9>c. 


(36) 


that is 


(q 2 j, q 2 k+i) = - (<? 2 fc+i, q 2 j) = rj8 jtk , {q 2j , q 2 k) = faj+i, 92 k+i) = 0. (37) 

These were first calculated in the paper [10] : 

q 2j (x)=x 2j , q 2j+ i(x) = x 2j+1 — 2jx 2j ~ 1 , rj = V2nT(2j + 1). (38) 

Given the choice of q /s described above, 

(a) the matrix A + B acquires a block-diagonal form 

A + B = r(8)J, r = diag(r 0 ,...,r n _i), J = ^ J ^ , (39) 

which leads to 


92n(z) 


Pf(r ( 8 ) J + (z - 1) A) 
Pf (r 0 J) 


(40) 
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(b) the matrix A is given by 


L 2j,2k 


= A 


2j+l,2k+l 


~ 0, A2j-l,2k — r ( j + k — - ) . 


(41) 


Notice that matrix elements of both r ® J and A labeled by a pair of indexes of 
the same parity vanish. Therefore, the 2 n x 2 n Pfaffians in the numerator and the 


denominator of (40) are reduced to n x n determinants: 

del ( f j 1 5j^ + (^ 1) A2j— l,2fc] l^j,k^n 


92n{Z) = 


det | kj [ ^jk 11 nj.kn : n 


(42) 


Finally, we apply the formula det(f/)/det(L 2 ) = det(V l UV x ) to perform divi¬ 
sion in (|42|). With the help of the explicit formulae (38) and (41) we get 


(*-i) ry + *-|) 

Ojk + 


\/ 2 i - l)r(2fc - 1) 


(43) 


- 1 ^j,k^n 


g 2 n{z) = det 
Lemma [l] is proved. 

A.2 Lemma 2. 

The proofs of Lemmas 2, 3 are based on the following integral representation for 
the matrix elements ([7]) of matrix M n : 


M n {j,k ) = 


dx 


x J 


x 


V2vr Jo x 5 / 2 V r ( 2 J - !) V r ( 2fc “ 1) 


, l^j,k< n, (44) 


which can be obtained by representing Y (j + k — 3/2) in ( Jt|) a s an integral. 
Take any v = (v\,V 2 , ■ ■ ., v n ) e M n \{0}. It follows from (|44[) that 


(v, M n v) = 


dx 


V2*Jo * 5/2 I U V r ( 2 i - 1) 


2 


VjX- 


. nr 3 


> 0 . 


(45) 


So, M n is positive definite by definition. 

Next, let us prove bound (J8| on the spectral radius of M n . Let Ai, A 2 ,..., X n >0 
be the eigenvalues of M n . Then 


/ a \ n 

= (A”„»)" < L AJ = (TrAC)" 


( 46 ) 


vfc= 1 
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It follows from the upper bound (13) of Lemma [3] that for any e > 0, there is N ( 
such that for any n > N e , 


TrIff™ ^ a/ —^ 7 1 nV —t £ — 1 — Ci + 6, 


where Ci e (0,1). Therefore, we can choose e small enough so that 


TrM™ ^ 1 - /i, 


where // e (0,1). Using this estimate in (46) for n > N e we get 

^max (j^) ^ (1 h) n U 1 . 


n 


Lemma [2] is proved for N = N e . 


(47) 


Remark. The spectral properties of M n seem quite interesting. For instance, 
in the large-n limit there is a concentration of eigenvalues near 1 such that the 
restriction of M n to the space spanned by the corresponding eigenvectors is close 
to an identity operator perturbed by an elliptic linear differential operator. Formal 
analysis of this perturbation suggests the asymptotic A mox (n) = 1 —/io'n . -1 + o(n _1 ) 
for suitable n o > 0 . 


A.3 Lemma 3. 


The integral representation (44) for the matrix elements of M n leads to the follow¬ 
ing integral representation for the trace of a power of M n : 


Tr M™ = 


dx i 


dx 9 


dx. 


™ _g—xi — X2—...—Xn 


Jo s/Znxi Jo \Z‘2nx2 Jo v 2 irx. 

COSh n _! (y/XmXx) COSh n _i (^/X 1 X 2 ) ■ . . COsll n _i ( y/X m -iX m ) , (48) 


where cosh„(x) = 2fc=o jkky. * s the degree-2n Taylor polynomial generated by the 
hyperbolic cosine. Performing the change of variables Xk = y\ hi (48) we can 
re-write the integral representation for Tras follows: 


TrM m = - 


7T 


m/2 


dye-^ y l 


cosh n _i( 2 / m j/i) cosh n _i(r/ir/ 2 ) • • • cosh n _i (y m -iyn 


(49) 


Here M+ = {(z/i, y 2 , ■ ■ ■, y m ) e | yk 3 s 0, k = 1, 2,..., m} is the first ‘quadrant’ 
of M m and dy is a shorthand notation for Lebesgue measure on M m . As the 
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integrand of (49) is symmetric with respect to reflection —> — y, for any i = 


1,2 ,, m, we can re-write Tras an integral over 

TrM™ = 7 / dye-^yl 


<=°shn-i(!M/i) c°sh„_i(j/iJ/ 2 ) - - -cosh„_i(y m _iy m ). 


(50) 


To prove Lemma 3 we will establish an upper and a lower bound on TrM " 1 and 
then compute the large-n limit of each of these bounds. 


A.3.1 An upper bound for Tr M™ 


A good starting point for the calculation is formula (50). For any ieK, cosh n _i (x) ^ 
cosh(x). Also, 

dz l-z~ 2n „„ 


cosh n _ 1 (x) = 


(51) 


2niz 1 — z~ 2 

where the integral is anti-clockwise around a circle of radius smaller than 1 centred 
at the origin in the complex plane. Replacing all but one cosh n _! with cosh we 
get: 


j dye s ^=i y * cosh n _i (y m yi) cosh(yiy 2 ) • • • cosh(y m _,y m ) 
^aia 2 ...a m _i [ dye~ E *- 1v * cosh n _! a iyiyi+i ? ( 52 ) 


where au, ct 2 , • • •, a m _i are independent identically distributed random variables 
which take values +1 with probability 1/2. Representing the remaining cosh n _i 


1 > 

m /2 

27T y 

) 

1 > 

m/2 

27T , 

) 


with the help of (51) and then computing resulting Gaussian integral over 
find 

m/2 


Tr AC < ( 2 


E, 


aia 2 ...am-i 


dz 1 — Z 


—2 n 


2niz 1 — z 2 




we 


(53) 


where 


/ 1 -si 
x 2 


0^(z) - del 


Ol | 
2 


0 

u 2 


02 

2 


0 

0 


\ 


0 

0 


0 

0 


0^3 
’ 2 


&m — 3 


■ £ 0 

2 u 


0 

0 


0 

0 


1 

— 2 


0 


Q m-2 

2 

l 

Q-m —1 


0 

0 


0 

Ot-m —1 


2 

i / 


(54) 
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The determinant can be calculated recursively in m, yielding D[ a \z) = 1 — z and 

m + 1 


D ( °\z) - -(m - 1)4- (z - A m ) [z + 


)<“>( 
for m ^ 2, 


(55) 


where A r 


n m—l 

k= 1 


Note that (55) implies that all principal minors of the 


matrix under the sign of the determinant in (54) are positive for z = 0. Therefore 


the matrix itself is positive definite for z = 0. By continuity, the real part of this 
matrix remains positive definite for z A 0 provided \z\ is small enough. Therefore, 
the real part of the quadratic form which determines the Gaussian integral in 


(52) is positive definite, which justifies the interchange of integrals leading to (53) 


provided the contour is taken to be a circle around the origin of a sufficiently small 
radius. 

Substituting (55) into (53) and changing the integration variable z —» A m z we 
find that the integrand no longer depends on ex’s. Averaging over ex’s becomes 
trivial and we get the following integral upper bound 


Trilf m < 


dz z 


—2 n 


27 rz z 2 — 1 Vl — z yj(m — 1 )z + rn + 1 


(56) 


The rest of the calculation is slightly different depending on whether m = 1 or 
m > 1. Here present the calculation for m > 1 only, the (simpler) case of m = 1 
can be treated along similar lines. We calculate the integral in the right hand side 
of (56) as follows. First we replace z~ 2n — 1 with z~ 2n in the integrand on the r.li.s. 


of (56), since this does not change the value of the integral as the omitted term 


is analytic inside of the contour of integration. Next we deform the contour away 
from the singularity at zero and out to infinity, leading to integrals around the 
other singularities of the (modified) integrand: a simple pole at z = —1, a branch 
cut singularity along the real line from 1 to +oo, and a branch cut singularity 
along the real line from — —oo. The contribution from the the integral 

over the large circle at infinity is zero. The contribution from the pole at z = — 1 
is easily evaluated as 1/4. Evaluating the integral around the branch from 1 to 
+oo it is convenient first integrate by parts, so that the singularity at z = 1 is 
integrable. The integrals along the two branch cuts lead to two real integrals 
whose asymptotics are controlled by the integrand (1 + y)~ 2n . Changing variable 
y —* y/Zn, and making some simple estimates on terms that do not affect the 
leading asymptotics, we are led to 


Tr Ml 


1 

< 4 + 


+ 


1 


n / 

vr m J o 
m + 1 


dy 

V^y 

m + 1 


yj'lnn t l \Jm - 1 V 2m 


1 + -*- 
2 n 

3/2 


—2 n 


(57) 


m — 1 
m + 1 


2n+l 


dy 

Vni 


fi + - 

V 2 n 


-2n+l 
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Both integrals in the above expression can be estimated using the following bound: 


Im = 


'o 


* JjL (i , ~ M 
V 7 u/ ' M - 


<1+ M’ 


(58) 


which follows by evaluating the integral, using the substitution t = (1 + ]fe) in 
terms of the beta function as 

Im 


- yf^-B (m - i 0 - Vmf(m - |)/r(M) 


and using bounds on the Gamma function. Using this in (57), the final result is 
Tr MU 


1 

55 4 + 


1 

55 4 + 


n 

77177 

n 

7rm 


n 





1 H-+ — a / — I 1 H- 


n 


8 V 7rn 




n 


m — 1 


m + 1 


2n— 8 



1 H-+ — a / — I 1 H- 


8 V 7rn 


n 


(59) 


which coincides with the claim (13) of Lemma |3j 

Dividing both sides of (59) by y/2n and taking the large n limit, we find that 


lim sup 


-=TrM- < 
V2 n 


277171 


(60) 


A.3.2 The limit lim n _»ao ^=Tr M™. 

The strategy is to derive an integral lower bound for Tr M™ and calculate the large 
n-limit of the bound. Our starting point is the relation (49) and the following 
estimate for the polynomial cosh n _!: 

Lemma 4. There exist two sequences (h n ) n ^ i, (S n ) n ^i a M. such that 
lim h n = \, lim S n = 2, 

n—>oo Z n—► oo 

e~ ny cosh n _i(ny) ^ h n l(y < S n ), for y ^ 0, n ^ 1. (61) 

Here t(y < S n ) is the indicator function of the set [0, S n ). 

In fact, as n —> oo, e~ ny cosh„_ 1 (?r?/) converges almost everywhere to (y < 2) 
for y ^ 0, but here we only need the lower bound. The proof of Lemma [4] is given 
in Section lA~4l 


Using the bound (61) in (49) we find that 

m/2 


/ o \ rn / z r 771 

TrilC n m ' 2 / dyV\Hyiy l+ i < S n )e-^^- y ^\ 

W Jr™ 7 = 7 


(62) 
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where y m +i '■= Vi- It is straightforward to verify that the domain of integration 


for the integral in (62) contains the hypercube (0, \/5V() m , 

Therefore, 


( 0 , vS) m <= {!/ e R’"\ytyt-n < s n , k - 1 , 2 ,... ,m}. 


n 1 (»< vs) < n HyiVi +1 < S n ), ye 


1=1 




(63) 


1=1 


Substituting (63) in (62) and changing the integration variables according to 


R = yi + V2 + + 

Zk = yk +1 -Vk, k = l,2,...,m- 1, 

we get the following lower bound: 


h m (‘l 

TrM™ A — \ - 


m/2 rm^fSn 

m/2 / di? 




'o 


(64) 


m \ir, 

/ d*!... dz m _ ie -f ejlt^+qslt 1 

where P m _i (/?,) is the intersection of the hypercube (0, \fSn) m and the hyperplane 

{y 8 -R+ | 2 /i + 2/2 + • • • + y m = R}- 


In the derivation of (64) we used the fact that the Jacobian of the transformation 


y —» ( R, z) is equal to 1/m. 

The large-n limit of the right hand side of (64) can be evaluated by arguing as 
in the Laplace method: 


liminf Tr M™ 

n—XX) y/2 n 

m/2 

. . . , n m/2 / dR 

n—>0 ° y2 nm \7ry J o ,/Rm-i 

n—>co V 2 V VT / 27T 


h m 

Y lim 


dz 1 ... dz^te-t Er-T^+QSLT^n 


n yy'm ^2 
7c 


= iimA/T h n{l 


m/2 


n 


(m- 1)/2 


= lim 


m/2 


n 


( m - 1)/2 


’ —00 

r 00 d\ ( r J 

L^AJJ ze 

v m 

2vr\ 2 /• co dA 




n 


2vr 


m \ 2 

e 


S n .„( 2 



= lim a -£K[- 
n-*oo V 2 V 7T 


m/2 


n(m -i)/2 1 M 2 


rt 


n J V 27rm V 27rm 


1 
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The crucial, albeit very standard, first step in the above derivation consists of 
verifying that extending the integration space for the z -integral from P m -i (/?.), 
when R e (0, 2), to M m_1 doesn’t change the large rt-limit. 

We conclude that 


lim inf 

n—>oo 


-=TrM n m ^ 
V 2n 


2irm' 


and in combination with (60) this gives 


lim 

n—>00 


V2 n 


2nm 


Statement (12) of Lemma [3] is proved. 


A.4 Lemma 4. 


Let (a n }“ =1 be an arbitrary sequence of positive real numbers which diverges as 
n —> oo slower than n 1//2 , that is lim a n = go, but lim a n n~ 1//2 = 0. We will show 

n —>oo n —>oo 

that there exists N 0 > 0 such that for any n > N 0 and x 5= 0 


: cosh n (nx) ^ \ \ - ^=a^e "” //4 ' N ) t(x ^ 2 — a n n 1 ^ 2 ). 

v47t ) 


(65) 


The statement of Lemm a |4| where cosh n (nx) is replaced by cosh n _i(mr), is easily 
deduced from equation (65). 

Our proof builds on the ideas of dedicated to the study of sections of ex¬ 
ponential series (Taylor polynomials generated by exp). Let e n be a section of 
exponential series defined by 


e n x = 




i — 1 7 

3 =0 J 

Consider also 

e^ +) (x) = e~ nx e n (nx) , e( l _) (x) = e~ nx e n (-nx). 

Then the function we are interested in can be written as 
fn (x) := e~ nx cosh n (nx) = ^ 

First we show that (x) > 0 for x ^ 0. One can check that 

L ( e ^e 2n (-lax)) = (: lux ) 2 " e 2 ” S 0. 
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( 66 ) 













and e 2nx e 2 n (— 2nx)\ x=0 = 1. So (x) 5= e 4na: > 0 for x ^ 0. The next step is 
to show that f n (x) is a decreasing function. However 


f'n ( X ) = ~ n An ( X 


which is negative by (66). 

The fact that f n (x) is decreasing and the positivity of (x) imply that for 
any non-negative x 


fn(x) ^ f n (x) t(x ^ 2 - a n n 1/2 ) 

2 * f n (2 - oi n rT l/2 ) l(x ^ 2 - a n n~ 1/2 ) 



Therefore, it remains to prove that 



(67) 


for all n > Nq, where No is chosen to satisfy a n n l ' 2 < 2 for all n > No. 

We start with a differential equation satisfied by . As it is easy to check, 


— e (+) (x) = _-_ 

dx n 1 ’ (n — 1)! 


' \n - 

nx) e 


( 68 ) 


So e ( n +) (x) is a decreasing function on 


Equation (|68|) has to be solved with a boundary condition lirrq,.^ (x) = 0, 

t n e~ nt dt. (69) 


which follows from the definition of The solution is 


ei +) (z) = 


n 


(n — 1 )! 


Let 


V27T n (n/e) r 


n\ 


By the Stirling approximation formula, (j) n = 1 + O (n 4 ) for n —> go and <f n < 1. 
Dehne 


t (t) = t — 1 — log t ^ 0, for t e 


In terms of <f n and r, expression (69) acquires the following form: 



(70) 


( 71 ) 
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The integral in the right hand side can be analysed using the Laplace method. It 
follows from the definition that 

pec 


1 = 4 +) ( 0 ) = n 



rW dt. 


0 


Therefore, (71) can be re-written as follows: 


3 (+) 


(x) = 1 - . 4> n e nT ^dt =: 1 — r n (x). 



ZL/l 


Let us estimate the remainder r n (x). Evidently, r n (x) ^ 0. An application 
of Taylor’s theorem with the Lagrange form of the remainder reveals that for 
0 < t ^ x ^ 1, 


T (t) ^ ^ (t — x) 2 + t' (x) (t — x) + T (x) . 


(72) 


1 — X 1 

Noticing that t' (x) =-and r" (x) = we can use the above bound on 

x x z 

r{t) to obtain the following upper bound on r n \ 


r n (x) ^ . <i> n e- nT W / e 2x 2 



n nil — x) 

it—x ) 2 +- d~ x ) 

X dt 


= (erfc (1 - x)j - erfc (2 - *))) 

s: ^xe-" T(l) erfex A /|, 

where erfc and erfex are complementary and scaled complementary error functions 
correspondingly. Finally applying the classical estimate erfex (x) ^ —— valid for 

xybr 

any x > 0 (see e.g. [I]) we obtain 


(x) < e -" r(a: ) 


'\f2mr 1 — x 
where we used that 0 n < 1. Therefore 


< 


x 


0 —nr{x) 


yphvK 1 — X 


Tin ( 1 - ~Y n ~ 


■ 1/2 ) < 

Using (72) for x = 1 and t = 1 — ^n -1 / 2 we obtain 

T 


I 1 — 1/2 

Tin ( 1 - Y n 




which leads to the desired bound (67) for . Lemma [4] is proved. 
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B On the numerical evaluation of p 2 n ,o- 

It is clear from the proof of Lemma [l] that the final form of the Pfaffian or determi- 
nantal expression for the probability that an n x n real Ginibre matrix has no real 
eigenvalues is strongly influenced by the choice of skew orthogonal polynomials 
used in the derivation. And even though the final exact result does not depend 
on the choice of the skew orthogonal polynomials, its numerical stability is highly 
sensitive to the choice. 

For example, the determinantal formula (|6]) is highly suitable for numerical 
evaluations since the condition number of the corresponding matrix / — M n grows 
at most linearly with n. Indeed, its largest eigenvalue is smaller than unity since 
M n is positive definite in virtue of the first part of Lemma [2} On the other hand, 
its smallest eigenvalue is separated from zero by an interval of length of order 
0(n _1 ) due to the result of Lemma [2] concerning the largest eigenvalue of M n . 

This should be contrasted to the determinantal formula derived in [13], The 
condition number of the matrix p appearing in this formula grows exponentially 
with n, forcing one to use high-precision numerics and leading to computation 
times growing exponentially with n. 
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